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First-principles calculations of the spin-orbit coupling in graphene with hydrogen adatoms in dense 
and dilute limits are presented. The chemisorbed hydrogen induces a giant local enhancement of 
spin-orbit coupling due to sp 3 hybridization which depends strongly on the local lattice distortion. 
Guided by the reduced symmetry and the local structure of the induced dipole moments we use 
group theory to propose realistic minimal Hamiltonians that reproduce the relevant spin-orbit effects 
for both single-side semihydrogenated graphene (graphone) and for a single hydrogen adatom in a 
large supercell. The principal linear spin-orbit band splittings are driven by the breaking of the 
local pseudospin inversion symmetry and the emergence of spin flips on the same sublattice. 



PACS numbers: 71.70.Ej, 73.22.Pr 

Spin-orbit coupling is central for a variety of spintron- 
ics phenomena [TJ [2] such as spin relaxation, spin trans- 
port, or topological quantum spin Hall effects. Itinerant 
electrons in graphene have weak spin-orbit coupling as 
they are formed from p z orbitals. The Dirac cones are 
separated by what is called intrinsic spin-orbit coupling 
of 2Ai = 24/ieV due to p z -d mixing [3-6 . This small 
value is desirable for long spin lifetimes, but experiments 
suggest [7HTQ] that spin relaxation is governed by much 
stronger spin-orbit, or perhaps magnetic [11 , interac- 
tions than the intrinsic one. Potential culprits are light 
adatoms [T2HT4] which are typically not important for 
momentum scattering but may be essential for spin-flip 
scattering. On the other hand, large spin-orbit coupling, 
when controlled, is desirable for engineering robust quan- 
tum topological phases in graphene covered with heavy 
adatoms [T5HT7] . 

Hydrogen is an ideal light adatom to study induced 
spin effects in graphene. Not only can it produce lo- 
cal magnetic moments [T8H2T] , as recently experimentally 
demonstrated [22., but it should also enhance graphene's 
spin orbit coupling (SOC), as proposed in Ref. [12j Unlike 
for heavy adatoms whose cores can directly contribute to 
SOC, the enhancement of SOC from hydrogen is solely 
due to sp 3 hybridization facilitated by local structural 
deformation. The SOC effects can thus directly probe 
sp 3 phenomena in graphene. 

The presence of both magnetic moment and large spin- 
orbit coupling makes the spin physics exciting, but also 
challenging to explain the spin relaxation experiments. 
To disentangle the two contributions, as well as to see 
what new phenomena they can lead to, it is important 
to consider them separately. Here we present a quantita- 
tive and qualitative study of SOC induced by hydrogen 
on graphene in two limits. One is the dense limit, rep- 
resented here by single-side semihydrogenated graphene 
(also called graphone) [23]. This structure is relatively 
simple and allows for a quantitative analysis of the sp 3 
hybridization effects on various spin-orbit parameters. 



Our results indeed show a giant enhancement of SOC, 
strongly dependent on the buckling deformations of this 
structure. We introduce a single-band and tight-binding 
Hamiltonians to describe the main SOC effects. In par- 
ticular, we show how pseudospin inversion asymmetry 
(PIA) introduces new terms (we call them PIA SOC), 
which couple the opposite spins on the same sublattice, in 
addition to Bychkov-Rashba nearest neighbor hoppings. 

We also quantify the local spin-orbit structure and pro- 
pose a minimal realistic SOC Hamiltonian in the dilute 
limit, represented by large supercells (starting with 5x5), 
intensively studied for orbital effects [24H27] . Based 
on first-principles calculations we demonstrate a giant 
local — and we identify the local impurity region from the 
dipole moments distribution — enhancement of SOC due 
to sp 3 hybridization and formulate a minimal realistic 
SOC hopping model. We believe that this is a bench- 
mark model to study spin relaxation, spin transport, but 
also weak (anti)localization [11, 28 phenomena in which 
spin- flip and spin-orbit scattering play an important role. 

Dense limit: single-side semihydrogenated graphene 
(SHG). The effects of sp 3 hybridization on SOC are 
studied using single-side semihydrogenated graphene [in- 
set in Fig.[TJa)] with different degrees of out-of-plane lat- 
tice distortion A of the hydrogenated carbon site which 
is on sublattice A. The C-H bond length g?h is 1.13 A, 
and we take the lattice constant to be the relaxed dis- 
tance a « 2.516 A between the nearest neighbors around 
the adatom site in the supercell calculations below; in 
graphene the lattice constant is 2.46 A. (The relaxed 
SHG structure of lattice constant 2.535 A would have 
A/a = 9.7% and the sp 3 tetrahedron 20.41%.) In exam- 
ples we choose A/a = 14% which corresponds to relaxed 
large supercell structures. 

To study SOC effects we restrict the computational 
basis to be spin unpolarized. The calculated electronic 
structure and the projected local density of states for 
SHG are shown in Fig.[]Jtop) for A/a = 14%. Compared 
to graphene, in which tt and 7r* bands without SOC touch 
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at K [29 J, the C-H bonding pulls them apart: the tt band, 
which at K comes mainly from sublattice A, is shifted to 
about 5 eV below the Fermi level (the GW approximation 
predicts a greater shift by about 2 eV [E]). The 7r* band, 
which comes from sublattice B, lies at the Fermi level. 
This band, which we consider for our SOC analysis, is 
narrow since the nearest-neighbor hopping is inhibited 
for p z electrons on B. The inset in Fig. [TJa) shows the 
7r* probability density at K that has the p z character on 
sublattice B. 
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FIG. 1. (Color online) Top: Calculated electronic band struc- 
ture of single-side semihydrogenated graphene. (a) Sublattice 
resolved band structure for the distortion A/a = 14%. The 
filled (red) circles correspond to sublattice A whereas the open 
(blue) to sublattice B. The circles radii correspond to the car- 
bon atom charge densities. The inset shows the structure and 
the probability density of the flat band at the K point, (b) Or- 
bital resolved local density of states. Bottom: Extracted spin- 
orbit coupling parameters for the 7r* band at K and T as 
functions of A/a: (c) Intrinsic spin-orbit coupling splitting 
2Ai at K. The inset shows the band splitting, (d) Adatoms- 
induced SOC splitting A§r at K. The inset shows the spin 
texture around K for the lower spin-orbit split band. The in- 
plane components are shown by the arrows while the z com- 
ponent by the color map. (e) Adatoms-induced SOC splitting 
Abr at r. The inset shows the band splitting around F with 
the identified Bychkov-Rashba energy Ebr = 0.87 /xeV for 
A/a = 14%. 



We now extract the SOC parameters for the states at 
K and T for the 7r* band which is at the Fermi level. 
Since this band is nondegenerate, the effective spin-orbit 



Hamiltonian can be expressed via the spin Pauli matri- 
ces s. The small group of K (T) is Cs (Cs v ). Up to 
terms linear in momentum, which is here measured from 
K and T, the SOC Hamiltonians compatible with those 
symmetries are 



H^ff — ^BR^x^y ~ k y s x ) + t\\s z 
^eff = ^BR^y^x — k x S y ) . 



(i) 

(2) 



Here r = 1 ( — 1) stands for K (K'), Ai is the adatom- 
modified intrinsic spin-orbit coupling and Abr is the 
adatom-induced (Bychkov-Rashba- like as in semiconduc- 
tor physics [30]) spin-orbit coupling. We will see that 
the latter comes from the space and pseudospin inver- 
sion asymmetry. Contrary to graphene, the BR SOC at 
K depends on the momentum magnitude. Higher-order 
terms in Wj§ and are presented in [3T] , 

Figure [T] (bottom) shows the extracted SOC param- 
eters as functions of A/a. The intrinsic SOC Ai is ob- 
tained from the splitting of the band at K, see inset in 
Fig. [ijc). Parameter Ag R is extracted by fitting the lin- 
ear dependence of the ratio of the spin expectation val- 
ues, (s x )/(s z )\ kx=0 = A B c R /c 2/ /Ai + 0(/c^), close to K. The 
trigonally warped spin texture around K is shown in the 
inset of Fig. j^d) and can be described by higher-order 
terms in T^jf [3T]. Finally, Ag R is obtained by fitting 
the spin splitting at T, see inset in Fig. [TJe). SOC is sig- 
nificantly enhanced in comparison to graphene. Directly 
comparable is the intrinsic SOC parameter whose value 
in graphene is 2Ai = 24 fieV [3]. 

The above single-band model can be obtained from 
a tight-binding (TB) Hamiltonian using the carbon p z 
and hydrogen s orbital basis. The Hamiltonian contains 
orbital and SOC parts, T~L = H rb + %so- We denote 
by c\ a = (a\ a ,b\ a ) and c ia = (a ia ,b ia ) the creation and 
annihilation operators for the p z orbitals on the sublat- 
tices (A,B), with spin a and lattice site i. Similarly, 
we define and h m(T for the hydrogen s orbitals on 
adatom sites m. For the orbital part H rb we take the 
TB model Hamiltonian introduced in Refs. |26] and [32j 
which assumes the nearest neighbor carbon-carbon (C- 
C) hopping t = 2.6 eV, direct carbon-hydrogen (C-H) 
hopping T, and the adatom on-site energy e^: 

n orh = e h ^ + T Yl h ™ c ia ~ tJ2 c ^ c j*' 

m (m,i) (ij) 

The angle brackets denote the nearest neighbors. Fitting 
the TB model to the first-principles band structure for 
A I a — 14% (distortion in the single- adatom limit, see 
below) we obtain = 3 eV and T = 6.5 eV. The values 
are reliable in the vicinity of K point where p z carbon 
orbitals dominate the projected local DOS. 

The SOC Hamiltonian can be derived by inspecting 
the reduction of the graphene point group symmetry 
Deh — which allows for the intrinsic spin-orbit coupling Ai 
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only — to the one corresponding to SHG. First, the C-H 
covalent bonds break the space inversion symmetry and 
the point group reduces to C$ v . This structure inversion 
asymmetry induces the Bychkov-Rashba-like term Abr. 
Second, the hydrogenated carbons on sublattice A cannot 
be interchanged with the non-hydrogenated carbons on 
sublattice B. This breaks the pseudospin inversion sym- 
metry and Cq v —> Czv- The effect of the latter reduction 
is twofold: (i) The intrinsic SOC depends on the sub- 
lattice: A^ and Af ; (ii) New SOC terms emerge due 
to the pseudospin inversion asymmetry, Ap IA and A^ IA , 
discussed below. 

As the hydrogen s orbitals do not directly contribute 
to SOC, we can express the SOC TB Hamiltonian in the 
p z basis. In the next-nearest-neighbor limit this Hamil- 
tonian has five real parameters and reads, 



K(K / ) we obtain, 



c ja , [A BR (axdy)J, 



E 



(id) 

J r 



ja' 



2A pia (s x D i. 



(4) 

The double angle bracket stands for the next nearest 
neighbors and label c denotes sublattice A or B. Fac- 
tors vij = 1 (—1) for clockwise (counterclockwise) hop- 
ping path j to i. The nearest-neighbor and next- 
nearest-neighbor Dij unit vectors point from j to i (in 
a flat lattice). The first term in Eq. Q is the stan- 
dard Bychkov-Rashba hopping as for graphene. The sec- 
ond term describes the sublattice resolved intrinsic SOC 
which couples same spins, and the PI A term which cou- 
ples opposite spins on the same sublattice. Hamiltonian 
Hso m Eq. Q applies to any hexagonal lattice system 
with Cs v point group symmetry, such as BN, or silicene 
in a transverse electric field. 

The single-band limit described by Eq. ([!]) can be ob- 
tained from the TB Hamiltonian H = H or b + ?4o by 
downfolding to sublattice B. This gives, 



A 



BR 

Ai 



-aA B 



PIA 



-Af 



aA B R ^ 

) A BR g h 



(5) 
(6) 



Both PIA and BR SOC hopping terms contribute to the 
effective band SOC parameters. This is the likely reason 
for the extracted nonmonotonic dependence of Ai and 
the decrease of Ag R as a function of A/a shown Fig. [I] 
(bottom), The TB model cannot be reliably used at Y 
as there other bands (orbitals) mix in, see Fig. [TJa) and 
Ref. El 

We also present an effective SOC Hamiltonian close 
to K. After transforming l~L so to the ordered Bloch ba- 
sis [ipA^(k)^ipAi(k)^B^(k)^Bi(k)] and linearizing near 



: A.BR(ra x s y 
- \ [Af +B cr z 



A 



A-B 



<T TS Z 



(7) 



\ [A^+ B a z + Ap IA B (j ] a{k x s y - k y s x ) . 



Here (ao, <x) and (sq? s) stand for the unit and Pauli ma- 
trices in the pseudospin and spin spaces, respectively. 
The momentum is measured form K(K') and parameters 
A A±B = A A ± A B . If z-inversion symmetry is restored, 
Abr, ApiA B > and A A_B vanish and one obtains the sil- 
icene limit [33] . 

Hydrogen on a supercell: single- adatom limit. The 
single-adatom limit is represented by a 5 x 5 supercell 
with a single hydrogen (2% coverage). We use a fully 
relaxed structure with A & 0.36 A (14% distortion) and 
the next nearest neighbor distance a — 2.516 A of the 
carbon atoms around the hydrogenated site Ch- 

Figure [2| a) shows the calculated spin-unpolarized elec- 
tronic band structure of our 5x5 supercell. The 
low energy spectrum contains three characteristic bands: 
the valence and conduction bands, and the mid-gap 
impurity-like band. These three bands can be nicely fit- 
ted by two parameters T = 7.5 eV, and = 0.16 eV en- 
tering the Hamiltonian H or b, Eq. as seen in Fig.[2|a). 
Our values differ from Refs. l26l and l32| the comparison 
is discussed in [31]. Larger supercells are also well de- 
scribed by these parameters, confirming that the 5x5 
one already describes the dilute limit, see [31]. The TB 
model with s and p z orbitals is also supported by the pro- 
jected p z local density of states per carbon atoms around 
the hydrogenated carbon as compared to the total den- 
sity of states, see Fig. [2^b). 

The calculated spin-orbit splittings are shown in Fig. 
|2|c). To explain them we propose a minimal realistic 
SOC model which is locally Cs v invariant in the impu- 
rity region with the C-H bond as the threefold axis of 
symmetry. We deduce the impurity region from the in- 
duced dipole moments, shown in Fig. |3ja). The main 
effects are confined up to the second nearest neighbors of 
the hydrogenated site Ch (in sublattice A), defining our 
impurity region. We use A\ (A a ) for the creation (anni- 
hilation) operators on Ch and G (B m jT ) on the three 
nearest neighbors. Otherwise the terminology follows the 
SHG case. The SOC Hamiltonian compatible with the 
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FIG. 2. (Color online) First-principles results and their tight-binding model fits for hydrogen adatom on 5 x 5 supercell. 
(a) Electronic band structure around the Fermi level. Dots are first-principles results and solid lines are tight-binding model 
fits. The band spin-orbit splitting around K is sketched for the valence band, indicating the out of plane z components and the 
rotation directions of the in-plane spin components, (b) Broadened total density of states per atom (gray) and p z projected 
local densities for atoms in the vicinity of the adatom: hydrogenated carbon atom Ch (green), its nearest neighbor C nn (dark 
blue), and the next-nearest neighbor C nn n (light blue), and the total density on hydrogen H (red). The projected densities of 
states are normalized to the corresponding number of atoms in the set. Conduction (c), impurity (d), and valence (e) band 
spin-orbit splittings along the high-symmetry lines; symbols as in (a). Tight-binding model least-square fits are performed 
within the shaded regions around K. (f) Spin expectation values around K for the spin-orbit split valence and conduction 
bands closer to the Fermi level. 



local symmetry is 

^so = g E C la C jv' 



Ai A 



E K <v 

«C H ,j» 



^3 



h.c. 



2i 



E A t B j<r* [ABR^Xdcn^J^+h, 



(8) 



<C H ,j) 



2i 



E< ^[APia^xD^J^- 



((id)) 

The first term is the graphene intrinsic SOC (2Ai = 
24/xeV). It couples all next-nearest neighbors pairs not 
containing (this is denoted by the primed summation 
symbol) Ch- The second term describes the adatom 
induced intrinsic spin-orbit coupling coupling Ai, which 
couples the same spins on the same sublattice. The third 
term, with Bychkov-Rashba hopping parameter Abr, de- 
scribes the induced nearest neighbor spin flips. Finally, 
the fourth term, with PIA parameter Apia, comes from 
the pseudospin inversion asymmetry. This term couples 
opposite spins of the next nearest neighbors. We remark 
that Csv symmetry allows more spin-orbit hopping terms 
in our impurity region. We considered them all but found 
only the three A's in Eq. ^ relevant to explain our ab- 
initio results, see the scheme in Fig. |3jb). 

Figures |2jc,d,e) show spin-orbit coupling induced band 




FIG. 3. (Color online) (a) First-principles calculations of elec- 
tric dipole moments induced by hydrogen adatoms on a 5 x 5 
supercell. Directions of the dipole moments are shown by ar- 
rows; the sphere radii correspond to the dipole magnitudes, 
(b) Hopping scheme of the tight-binding model showing the 
relevant orbital and spin-orbit coupling parameters. 



splittings along high symmetry lines. The multiband 
least-square fit around K point gives the following val- 
ues for the SOC parameters: Aj = —0.21 meV, which is 
about 17 times larger than that of graphene Ai; Abr = 
0.33 meV, more than 60 times the value in graphene 
where Abr = 5 /ieV in a representative transverse elec- 



tric field of 1 V/nm [3]; A^ IA = —0.77 meV, which has 
no counterpart in fiat graphene. The signs of the above 
parameters have been determined from the spin expec- 
tation values around the K point, shown in Fig. [2^f ) . 
The spin texture is governed mainly by PIA SOC. Those 
parameters also fit larger supercells [31] and we propose 
them, together with Hamiltonian Eq. (JsJ, to describe the 
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single adatom limit important for investigating spin-flip 
and spin-orbit scattering in graphene. 

In conclusion, we investigated spin-orbit coupling in- 
duced by hydrogen, representing light adatoms, on 
graphene in dense and dilute limits. We introduced real- 
istic model spin-orbit Hamiltonians and provided quan- 
titative values for their parameters that can be used to 
study spin relaxation, spin transport, and mesoscopic 
transport in graphene with adatoms or in similar two- 
dimensional structures of the same symmetry. 

We thank T. O. Wehling for helpful discussions. This 
work was supported by the DFG SFB 689, SPP 1285, 
and GRK 1570. 
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Computational method 

The electronic structure calculations are based on den- 
sity functional theory [1]. Generalized gradient approxi- 
mation [2] has been applied to the exchange-correlation 
potential and used in the full-potential linearized aug- 
mented plane wave (FLAPW) method as implemented in 
the FLEUR code [3]. The self-consistent calculations with 
spin-orbit coupling (SOC) in the film mode for semihy- 
drogenated graphene were carried out with cut-off pa- 
rameter k max of 6.7 bohr -1 and 144 fc-points in the ir- 
reducible wedge of the Brillouin zone. For the super- 
cell calculations we used k max — 4.7 bohr -1 and 64 k- 
points in the irreducible wedge of the Brillouin zone. The 
muffin-tin radii for carbon of 1.32 bohr and for hydrogen 
0.81 bohr were taken. 

The atomic positions in the supercell calculations have 
been relaxed using the quasi-newton algorithm based on 
the trust radius procedure implemented within the plane 
wave pseudopotential code Quantum Espresso [4]. For 
carbon and hydrogen we have used ultra-soft pseudopo- 
tentials [5, 6] with PBE exchange-correlation functional 
[2] with kinetic energies cut-offs of 30 Ry for the wave 
function and 300 Ry for the density. We used vacuum of 
15 A to separate the hydrogenated graphene planes. In 
Table I we list the relevant structural distances around 
the adatom site for different supercell sizes n. Covalent 
C-H bond length dn is close to 1.13 A. Next-nearest 
distance a between the three closest carbon atoms to 
the hydrogenated carbon, a tetrahedral edge length, is 
2.516 A. The hydrogen bond distorts the graphene plane 
by pulling the hydrogenated carbon atom out of plane. 
Parameter A corresponds to the vertical shift of the hy- 
drogenated carbon atom with respect to its three nearest- 
neighbors. The ratio A/a is about 14% in all stud- 
ied supercells. We note that for an ideal tetrahedron 
(methane) the ratio A/a — 20.41%. The binding energy, 
£t>ind = Eq — Eq 11 — Eq, of the hydrogen on graphene, 
calculated as the difference between ground state energy 
Eq of the hydrogenated graphene and ground state en- 
ergies of graphene without hydrogen Eq 11 and isolated 
hydrogen atom Eq is about —2.9 eV. 

The fine structure of relaxed atomic positions, espe- 
cially the local distortion in terms of the vertical shift 
A of the hydrogenated carbon atom, see Fig. 1, signifi- 
cantly affects the impurity band width and consequently 
the value of the on-site impurity energy parameter en- 



n x n 


du [A] 


a [A] 


A [A] 


A/a 


£bind [eV] 


4x4 


1.1286 


2.514 


0.359 


14.264% 


-2.83 


5x5 


1.1279 


2.516 


0.363 


14.415% 


-2.88 


6x6 


1.1283 


2.516 


0.362 


14.399% 


-2.87 


7x7 


1.1277 


2.517 


0.360 


14.297% 


-2.92 


8x8 


1.1281 


2.516 


0.361 


14.362% 


-2.94 


9x9 


1.1274 


2.517 


0.353 


14.012% 


-2.86 


10 x 10 


1.1283 


2.516 


0.360 


14.303% 


-2.91 



TABLE I. Structural parameters for supercell calculations ob- 
tained using Quantum Espresso code. Covalent bond carbon- 
hydrogen length dn, distance between the carbon atoms 
around the hydrogenated carbon a, vertical distance between 
hydrogenated carbon and its nearest neighbors A, and bind- 
ing energy -Ebind are listed. 

tering the tight-binding model description. The more is 
the hydrogenated carbon pulled out of plane, the wider is 
the impurity band and the larger is needed to describe 
the impurity band. 




FIG. 1. Local distortion due to chemisorbed hydrogen on 
graphene gives the covalent bond of length du between hy- 
drogen and hydrogenated carbon, increased distance between 
the next-nearest neighbors around the hydrogenated site a, 
and pulling the hydrogenated carbon out of plane by A with 
respect to the three nearest neighbors of the hydrogenated 
site. 

Performing least-square fits to the ab-initio data 
around K we found T = 7.5 eV and = 0.16 eV. Solid 
lines in Fig. 3(a) in the paper show the model fit to the 
5x5 supercell ab-initio data. In Fig. 2 we demonstrate 
the robustness of the model parameters by comparing 
with ab-initio data of larger supercells. The larger is the 
supercell the more narrower are the pseudogaps and the 
more flat is the impurity band. In the limit of an infinite 
supercell the band structure should contain the hydrogen 
s energy level and the graphene band structure with the 
Dirac cone. We note that previous studies [7, 8] present 
different values for the C-H hopping T and the on-site 
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FIG. 2. Calculated electronic band structures of hydrogenated graphene with different supercell sizes (hydrogen concentrations). 
The first-principles data are shown by (gray) dots whereas the tight-binding model fits by (blue) solid lines. The tight-binding 
model band structure with parameters from Ref. [7] T = 5.2 eV, £h = 0.26 eV and Ref. [8] T = 5.72 eV, £h — 1.716 eV are 
shown with (red) dashed lines and (blue) dashed-dotted lines, respectively, (a) 7x7; (b) 9 x 9; (c) 10 x 10 supercells. The 
bottom figures are detailed views of the impurity band for each supercell. 



energy e^] the comparison is shown in Fig. 2. 

To calculate the dipole moment of an atom A at atomic 
center f^, we used pa (ta) — J Ra ^(r y )(r / — f^d 3 ?"', where 
Ra is the atom muffin-tin radius, and g(f) is the va- 
lence electronic charge. The integral is performed using 
a weighted sum for an arbitrary polynomial accuracy [9] . 
For this purpose we have implemented the new module 
r3dgl into FLEUR. 

In Fig. 3 we show the dipole moment magnitude as 
a function of the in-plane distance from the adatom for 
a 5 x 5 supercell. There are significant dipole moments 
on the nearest-neighbor carbon atoms around the hy- 
drogenated carbon. The dipole moments diminish at a 
distance of about 4 A. In the inset to Fig. 3 we plot 
the in-plane components of the dipole moment directions 
(unit vectors) . This figure quantifies the dipole-moments 
distribution picture in Fig. 2 of the paper. 
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FIG. 3. Calculated spatial distribution of dipole moments of 
a 5 x 5 supercell with one hydrogen adatom. The magnitude 
of the dipole moment is shown as a function of the in-plane 
distance from the hydrogenated carbon. Inset shows the plane 
projection of the dipole moment directions. 



Single band spin-orbit Hamiltonian: symmetry 
approach 

Here we focus on the single-band spin-orbit effective 
Hamiltonian valid near K(K') point. The guiding princi- 
ple for constructing such a Hamiltonian is the invariance 
principle with respect to the small group of K(K'). 

The group symmetry for both the single top positioned 
impurity and the single-side semihydrogenated graphene 
(SHG) is Cs v . The threefold axis (rotations 7Zo — Id 
and 7Z±2tt/3 about the z-axis) for the single impurity 



graphene is the carbon-hydrogen bond and three reflec- 
tion planes (a v and a v o7^±2tt/3) are defined by the three- 
fold axis and the lines connecting the hydrogenated car- 
bon with its three nearest neighbors. The arrangement 
is slightly different in SHG, which has the threefold axis 
in the center of the primitive hexagonal cell, and the re- 
flection planes are defined by the threefold axis and the 
lines connecting the opposite vertices of the hexagon, see 
Fig. 4. 

The Brillouin zone remains hexagonal. The rotational 
subgroup Cs C C^v leaves K and K 7 invariant. The min- 
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FIG. 4. Coordinate system and symmetry operations of the 
point group for single-side semihydrogenated graphene. The 
hydrogen atoms are not shown. 

imal rotationally invariant 2x2 Hamiltonian reads 

Uf s = [a + bA: + A:_+c(4 + /c 3 _) + ...]^ + ^ 

+ [a k- + b k\ + c k + k 2 _ H ] s + + H.c. . 

Here k x and k y are momenta measured from K or K 7 , 
k± — k x ± iky, and s± = \{s x ± zs y ) are the spin raising 
and lowering operators; a, b, c, . . . are real and a, b, c, . • • 
complex parameters. Applying the time-reversal opera- 
tor T to we get an effective Hamiltonian Tif^ valid 
near K 7 , but K 7 can be mapped back to K using a vertical 
reflection o v £ C3^, i. e., apart from the rotations also 
the operator a v o T belongs to (the spinor extension of) 
the small group of K(K'). The small group invariance 
constraint, 

(<7„oT)Kff =K« K°T), (2) 

restricts a, b, c, . • • to be purely imaginary, i. e., a = ia, 
b = i/3, c = ij, . . . The expansion in momentum can be 
increased to an arbitrary order and can be rewritten in 
the form 

H T e i=n(k x ,k y ,r)-s, (3) 

where the valley index r = 1 (— 1) for K (K 7 ) points. The 
spin-orbit field n(k x , k y: r) up to the third order reads 

n x = ak y - rf3{2k x k y ) + -f(k 2 x k y + fc;}) , 

n y = -ak x + Tp(k 2 y - hi) - j(k x + k x k 2 y ) , (4) 

n z = ra + rb(k 2 x + k 2 y ) + c(2/^ - 6k x k 2 y ) . 

From the effective Hamiltonian Wj^ one can deduce 
basic properties of the underlying spin-orbit field tex- 
ture in the vicinity of K(K') point. For example, close 
to K(K') the spin-orbit field is circulating, being of the 
semiconductor Bychkov-Rashba-type but with the valley 
dependent spin ^-component, 

a(k y s x - k x s y ) + r &s z . (5) 



The above formula is used in the paper (a = Ag R , a = Ai) 
with SOC parameters extracted from the spin-orbit field 
around the K point, see Fig. 1 in paper. Away from 
K(K ; ) nonlinear terms begin to play role, Eqs. (4), and 
the spin texture becomes trigonally warped, see Fig. 5. 



(a) (b) 




FIG. 5. Calculated spin-orbit effects along a contour around 
the K point for the lower energy spin-split ir* band of SHG. 
The out of plane spin component for the spin-split band un- 
der consideration points opposite to z direction, (a) Spin 
expectation values along a contour whose radius is 1% from 
the Brillouin zone width, (b) Polar plot of the in-plane 
spin amplitude, a/ (s x ) 2 + i s y) 2 - The solid line is a model 
fit to the ab-initio data including the quadratic terms in 
Eqs. (4); the fitted SOC parameters are a — 0.102 meVnm 
and (3 = -0.119 meVnm 2 . 

Similarly, other high symmetry points in the Brillouin 
zone can be investigated. For example, in the vicinity of 
T the SOC Hamiltonian reads 

^eff ~ ^br(^^ — k x s y ) . (6) 

As discussed in the paper, the single band model can 
be microscopically justified by a tight-binding model with 
the basis of p z carbon and s hydrogen orbitals. In Fig. 6 
we show the spatial distribution of the probability density 
for a SHG at K and T points for the tt* band close to the 
Fermi level. From the spatial distribution one recognizes 
the dominant p z character of the states which supports 
our tight-binding model built on the p z basis. 




FIG. 6. Probability density plots in single-side semihydro- 
genated graphene for band 7r*, close to the Fermi level, for 
(a) K point and (b) V point. 
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